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Summary: The thcrmodynamical phenomenological theory of magnetically active superconducting 
materials and magnetic-superconducting heterostructures is presented. The materials may exhibit arbitrarily 
strong anisotropy, parametric or structural. Exact anisotropic similarity transformations are found for both 
the continuum and layered hybrid systems. 
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1 Introduction 

A formulation of thermodynamical theory of magnetically active, anisotropic materials ad- 
mitting coexistence of the superconductivity and magnetic order has been proposed in 
The theory is based on the Ginzburg-Landau approach extended to multi-component order 
parameters and the states of thermodynamic quasi-equilibrium far from the superconduct- 
ing phase transition. The field equations are derived under assumption of the U{1) gauge 
invariance. The nonlinearity with respect to the order parameter, magnetic field, and other 
physical quantities is assumed reasonably general, without any particular restrictions of its 
functional form. 

In the present paper we give an improved version of that theory complemented with the 
exact anisotropic similarity transformations for layered hybrid systems [23]. 

2 Background 

A material medium is considered magnetically active if its magnetic properties differ (sub- 
stantially) from those of the physical vacuum. Although the superconductivity itself is 
a magnetic phenomenon, the superconducting materials need not be magnetically active. 
This circumstance justifies disregarding the magetic activity of a superconducting material 
in heuristic or didactic considerations, but not in a more complete theory. Due to orbital 
polarizability and/or spin degrees of freedom, the effects of magnetic activity of superconduc- 
ing materials are practically always present, and in many instances significant. In particular 
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it concerns the heterostructures which involve interactions between superconducting and 
ferromagnetic components [6l [18]. One may also notice that, on the background of some 
controversial views [H |H] concerning coexistence of superconductivity and ferromagnetism, 
materials which combine these phenomena were found experimentally [21 [HI [H [15] • More- 
over, it seems that, in accordance with earlier theoretical predictions, the superconductivity 
of such materials is magnetically mediated [12] . 

In order to describe the superconducting state of condensed matter, the GlNZBURG- 
Landau theory [H] introduces a phenomenological quantity, called the order parameter, rep- 
resented by a complex scalar field '^(a:) with properties quite analogous to the microscopic 
quantum-mechanical single-particle wave function. Contrary, however, to the microscopic 
wave function, the GL order parameter describes the macroscopic condensate of the super- 
current carriers and enters as an independent variable into the thermodynamical free energy 
functional. This approach turned out to be extremely succesful. 

In the original GL paper the functional was given a simplest form suitable for description 
of the states of a superconducting material near the phase transition. In this form, the 
theory was confirmed afterwards by GORKOV [10] through rederivation from the microscopic 
theory. Its range of applicability turned out, however, severely restricted by the condition of 
proximity to the transition point. 

To cope with broader ranges of real situations, the GL theory has been extended in 
various directions, including multicomponent order parameters, more elaborate forms of 
the free energy functional, time-dependent phenomena, and heterostructures. Any such 
modification, as long as it respects the idea of a (possibly generalized) superconducting 
order parameter as a thermodynamical variable, is referred to as belonging to the framework 
of extended Ginzburg-Landau phenomenology [211 1221 [23] • 

The theory under consideration admits arbitrary (in a reasonable sense) functional forms 
of the GL thermodynamical potential. We intend to include in this way also the thermody- 
namical states which are far from the superconducting phase transition. We admit multiple 
order parameters, anisotropy, and material inhomogeneities, but, for simplicity, we retain 
the assumptions of stationarity, thermodynamical quasi-equilibrium, and charge neutrality. 

The assumption of charge neutrality refers to the resultant electric charge density of all 
charged components, such as the superconducting carriers (Cooper pairs), normal electrons, 
and the ionic substrate. It restricts our considerations to the case of bulk superconductivity 
in good (quasi-metallic) superconductors. To include "bad" superconductors and interface 
effects one would need a further extension of the model, which is beyond the scope of the 
present paper and will be given elsewhere. 

We shall focus our attention on superconducting materials which can be (a) strongly 
anisotropic, (b) inhomogeneous at the macro- and micro-scale, and (c) magnetically active. 
Below we shall comment briefly on the above marked specific assumptions. 

(a) The idea of isotropic superconductor, due to its simplicity and heuristic power, at- 
tracts much attention as a tool in theoretical investigations. It suffices here to note examples 
of most outstanding and successful classic theories of superconductivity such as Ginzburg- 
Landau, Bardeen-Cooper-Schriefer, and Abrikosov's ones, each of them based on an appro- 
priate isotropic model. On the other hand, the real superconducting materials exhibit, as a 
rule, some anisotropy of physical properties. In particular, the anisotropy of high cuprate 
oxides is rather high [SI [7]: e.g. the ratio of penetration depths Xc/^ab for YBa2Cu307_2^ 
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and Bi2Sr2CaCu208+j/ lies in the range 5-^7 and 50-=-200, respectively; it corresponds to 
estimated values of the ratio of flux line tensions Uc/Uab ranging from 20 to 40000. 

(b) To describe heterogeneous superconducting materials or devices, the theory must 
admit material inhomogeneities of various - macroscopic, microscopic, and intermediate - 
scales of length. The inhomogeneities can be smooth or abrupt, technological or natural. 
Many high Tc superconducting materials exhibit a distinct layered structure at the atomistic 
level. The spacings between atomic layers, such as Cu-0 planes in cuprate-oxide supercon- 
ductors, can be relatively large as compared with the in-plane distances what, besides the 
structural anisoptropy, can result in a natural micro-scale inhomogeneity of those materials 
in the out-of-plane direction. The Lawrence-Doniach model [13] and its generalization 
[22| [23] can be considered an extreme case of such a natural micro-inhomogeneity. On the 
other hand, due to the recent progress in nanotechnology it is possible to fabricate diverse 
nano-scale technological heterostructures. 

(c) Apart from magnetic activity originating from physical factors, we found it also ex- 
pedient to introduce its equivalent at the level of theoretical modelling. We refer here to the 
anisotropic similarity transformations - a useful tool in the study of anisotropic supercon- 
ductors, discussed in Section ID 

3 Notation 

Whenever appropriate, we use implicit tensor notation with supressed indices; the indices are 
exposed only when otherwise the expressions would be equivocal or illegible. For instance, a 
point in the physical space (as well as the corresponding material point when the distinction 
makes no difference) is denoted simply by x, which is to be developed to (x^, x^, x^) in typical 
3D situation. 

Basically, the standard notation from recent literature on phenomenology of supercon- 
ductivity [20] is employed throughout the paper. The effective mass m* of the supercurrent 
carriers is tensorial; in the case of transversally isotropic superconducting materials (e.g. lay- 
ered oxides, in good approximation) the relevant mass tensor components are ruab (in-plane) 
and rric (out-of-plane). 

The complex conjugate of the order parameter ip is denoted by ip. The symbols d and V 
are employed in a somewhat non-standard way. The partial derivative operator is denoted 
by d, while the symbol V is reserved for the electromagnetic f/(l)-covariant derivative. Both 
operators d and V act from the left. We found it expedient to introduce also the conjugate 
operators d and V which act from the right. 

Due to a bug in the routine translating from TeX to PDF, in the final printed form of 
the reference [1] a repeated typographic error occured. The error consisted in the fusion of 
two consecutive short overbars into a longer one. It resulted in erroneous replacement of the 
expressions d with ipd and ipV with ipV. The present text contains terms counteracing 
those faults. 
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4 The reference model 



To fix the point of departure for construction of extended models we start from the reference 
model which, excepting the details of notation, almost exactly coincides with the original 
GL model [9]. It is also convenient to develop the notation adapted in this paper. 

1. The superconducting state of the superconductor is described with the aid of the order 
parameter represented by a single- component complex scalar function ip = ip{x). 

The gauge transformation for this parameter is defined as 

A + dA, 

^ ^ e^V, (1) 
ijj — )■ ipe^i^^, 

where e* denotes the electric charge of the supercurrent carriers. The corresponding covariant 
derivatives V and V are given by 

v^a-^A, (2) 



and 



V = a+^A (3) 
nc 

2. The free energy of the superconductor is given by the functional 

F[T-ij,A]= fd'xFs (4) 



with an appropriate free energy density Fs of the superconducting state. 
3. The free energy density Fs equals 

Fs = F^ + a\ij\'+^\i:\' VV^, (5) 

where F^ refers to the normal state and m* represents the effective mass of the supercurrent 
carriers. The coefficients a and /3 are given as 

a = a{T — Tc), a = const, (3 = const. (6) 

The normal state free energy density F^- does not depend on the order parameter. 

The expression ([5]) ensures the already required invariance with respect to the gauge 
transformations ([1]). 



5 Local first order functional 

The reference model sketched in the previous section relies upon many intentional simplifi- 
cations, including isotropy, homogeneity, and proximity to the N-S transition point. Some of 
these drawbacks are direct consequences of the particular form of the free energy density Fs 
and can be eliminated simply by admitting, in place of ([5]), more general functionals defined 
on the same space of the field variables. In this way one can successfully take into account 
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many realistic features of superconductors. As examples one can mention such properties 
as high anisotropy or more subtle thermodynamical behaviour of high superconductors 



The range of applicability of such models is, however, severely restricted by the singlet 
nature of the order parameter. From the microscopic point of view a singlet order parameter 
is justified, first of all, for superconductors with dominating isotropic s-pairing. Also other 
types of pairing can be treated in this way, provided that a single pairing mode, corresponding 
to a one- dimensional representation of the internal symmetry group, is strongly dominating. 

Many superconding materials require, however, multiple order parameters. One can men- 
tion here the heavy fermion unconventional superconductors [17] which exhibit anisotropic 
spin triplet p-pairing, possibly mixed with spin singlet s-pairing. Another typical example is 
delivered by the coexistence of spin singlet s- and c?-pairings in orthorombic high Tc materials 



In the present section we shall briefly discuss an extension of the GL phenomenology 
which avoids the above mentioned drawbacks. 

1. First of all we modify the assumption stated in item 1 of Section 1. In place of the 
singlet order parameter, represented in the reference model by a complex- valued scalar field, 
a more general multi-component order parameter is introduced. For brevity, this order pa- 
rameter will also be denoted by the single symbol ip = ip{x). In geometrical interpretation 
the admissible values of ip will be considered as organized into a differential manifold. When- 
ever a particular component will be referred to, an indexed symbol like = ip^{x) will be 
used and interpreted as an appropriate - real or complex - coordinate in the corresponding 
local map. The index A will take integer values from 1 to A^, with arbitrary fixed A^. The 
real dimension of the order parameter manifold may, due to the presence of complex- valued 
coordinates, differ from A^. 

2. From the physical point of view the multiplet order parameter represents generally 
a mixture of charged fields. In contrast to the reference model, such a mixture can not be 
characterized by a single scalar charge. To characterize the coupling of the multiplet order 
parameter with the electromagnetic field we must specify an adequate gauge transforma- 
tion. To that end, in place of the scalar charge e* we introduce a real symmetric matrix 



such that the gauge transformation is defined by the corresponding matrix reinterpretation 
of the old formula ([1]). In this way also the gauge covariant derivatives ([2]) and ([3]) retain 
their graphic form with reinterpreted meaning. 

In general, the charge matrix e* need not be diagonal. However, being real and sym- 
metric, it can be diagonalized by an appropriate orthogonal transformation of the order 
parameters ip"^. The eigenvalue 0, if present, corresponds to electrically neutral subspace of 
order parameters which carries no supercurrent; in spite of this property such order param- 
eters can seriously affect the thermodynamical properties of a superconducting system. The 
remaining eigenvalues of e* determine the spectrum of charges of supercurrent carriers. 

In most cases the supercurrent carriers are experimentally identified as electrically equiv- 
alent to pairs of electrons or holes. Then the spectrum of the matrix e* simplifies to at most 
3 eigenvalues: ±2e and 0. 



m 




(7) 
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3. We accept the assumption formulated in item 2 of Section |U so that the free energy 
functional retains the local form We, however, modify item 3.: instead of the particular 
form (E]) we assume 

Fs = F[,lj, A] = F{T, X, ^, d^, ^ B, A, dA). (8) 

The explicit dependence on x is intended to reflect natural or artificial inhomogeneities 
of the materials, such as inhomogeneous doping, structure defects, etc., including possible 
heterostructural properties of the system. Implicit dependence on some external control 
parameters (such as mechanical deformation) is admitted; the implicit parameters can also 
depend on the position x. 

Apart from reasonable mathematical properties and the gauge invariance, no particular 
restrictions will, in general, be imposed upon the function ([H]). The symbol e* stands here 
for the matrix ([7]). 

4. The requirement of gauge invariance results in the following relations 

OF - dF OF - - OF 

dF i dF - dF 

and 

OF dF 

d{d,Ak) ^ d{duA,) ~ ^ ^ 

to be satisfied identically by the function ([H]). 

The conditions (fTOj) and ( ITTi) can be integrated in a straightforward way. As a result, 
one obtains the functional dependence 

Fs = A] = F{T, X, ij, ^, ^ V, B) (12) 

with the vector potential A entering indirectly through the gauge covariant derivatives V, V 
and the (gauge invariant) magnetic induction field B = curl A. The gauge invariance condi- 
tions for the function ( fT2l) reduce to 

dF ^ ^ - dF dF 7^ * dF 

5. Taking into account the functional dependence ( lT2l) . the functional derivatives of Fs 
with respect to the field variations Sip, Sip, and SA can be represented in the form 

SF _dF dF 

Sip ~ dip d{\/ip)' ^ ^ 

SF _dF 9F 

Sip ~ dip d{ipV) ' ^ ^ 

and 

SF i dF ^ , , dF 

ja = Yc^"^' dm - wff + ^- ^ ^ 
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6. In consequence, one obtains the generalized Ginzburg-Landau equation 



dip d{i) V) 

and the magnetic field equation 



OF OF - 

V = 0, (17) 



with the intensity of magnetic field 
and the supercurrent density 



curl H = — j 

c 



dF 

H = 4n- (19) 



?■ - BF dF 

7. Taken literally, the above equations are derived for sufficiently smooth fields. As a 
rule, such degree of differetiability can not be guaranteed at abrupt interfaces, in particular 
at the heterostructural ones. However, the integral form of those equations is valid under 
much weaker differentiability conditions, sufficient for the derivation of boundary conditions. 
The standard procedure applied to the equation (fHI) leads to the boundary conditions, which 
in index notation take the form of continuity equation across the interface 

dF' dF'' ^ 

where n'^ represents the normal vector of the surface demarcating the regions 1 and II. In 
consequence of this equation the normal component of the supercurrent fl20l) is continuous. 



6 Layered hybrid model 

The free energy functional for higher-grade hybrid model [23j has the following form 

T = To + J's + I giB)£x. (22) 

The last term represents the contribution of magnetic field to the free energy. In the standard 
particular case of magnetically inactive material the function g(B) reduces to B^/Svr. The 
term J^q describes the normal state, while J-g the superconducting one. The supercoducting 
term is composed of two parts: 

•F5 = J'p + ^j, (23) 

where the part 

J^p = f dxdyFn (24) 

describes the contribution of atomic planes, while J-'j corresponds to interplanar Josephson's 
bonds. Fn has the Ginzburg-Landau form ([5]), with 2D operator V and the coefficients 
dependent, in general, on the plane n. 
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The form of the functional J-p already ensures its invariance with respect to the gauge 
transformation ([1]), understood this time in the 2D sense. 

The Josephson term in the free energy functional for the hybrid model of grade K has 
the form 

J^j = sJ2J2 / ^^(^y ^in- (25) 

The quantity Egn denotes the energy of J-link between n-th and (n+q)-th planes and is given 
by 

^gn = 2{Cg(IV^n|' + l^n+gH " 7g (V^nV^n+gC-^^''" + C.C.)}, (26) 

with the coupling parameters (g and 7^ vanishing for q > K. The exponent Pgn is defined 
by the formula 

g* r{n+g)s 

Pqn = T- / ^zdz. (27) 

lie J ns 

7 Structured order parameters, nonlocal and higher 
order models 

In the extension of GL phenomenology discussed in Section [5] we restricted our attention 
to superconductors which can be described by local first order functionals with multiplet 
order parameters. Such a restriction greatly simplifies the theoretical considerations. On 
the other hand, it eliminates some physical phenomena which are known to be present in 
superconductors. It concerns, in particular, the nonlocal effects and nonparabolic dispersion 
curves. Such phenomena can, however, be taken into account by admitting yet more general 
- nonlocal and/or higher order - free energy functionals. 

The microscopic theories of superconducting materials of complex atomic and electronic 
structure indicate various possible mechanisms of superconductivity [2]. To cope with the 
diversity of those mechanisms, the phenomenological description can be further generalized 
by introducing the structured order parameter. Instead of the customary order parameter 
represented by a (possibly multicomponent) function ip{x), we introduce a generalized one, 
represented by a function of '?/'(x, ^), where the new variable ^ runs over a certain internal 
space S. The space S serves as a mathematical stage on which the internal structure of 
superconducting carriers (such as Cooper pairs) appears. 

The space S depends on the material under consideration and, parametrically, on the 
temperature, doping, and other control parameters such as external mechanical deformation. 
In the case of uniform material samples the space S is the same for all material points; in the 
case of heterostructures and inhomogeneous material samples it can become x-dependent. 

The concept of structured order parameter offers tighter bound between the phenomenol- 
ogy and the microscopic theories. Apart from more accurate quantitative description it al- 
lows also to study possible topological transitions relalated to changes in the deep topological 
structure of the space S. The detailed discussion of those effects is, however, beyond the 
scope of the present paper. 
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8 Anisotropic similarity and scaling 




Assume now that the order parameter space is endowed with a hnear structure and consider 
a hnear transformation x ^ x' ,ip ^ ip' , A ^ A' of the following form 



(28) 



where Q, R, and Q are linear matrices; the matrices are real except R which can be, totally 
or partially, complex. The transformation rule 

d' = OQ-' (29) 

is an immediate consequence of (1251) . 

1. The transformation (128|) in its general form violates the U{1) gauge structure. If 
we request that, under (128|) . the depart gauge covariant derivative be transformed into the 
target gauge covariant derivative, V'?/' — )• Vip', the transformation (!28|) must be adequately 
resticted. It follows that the relation Q = must hold. Such a transformation will be 
called a (gauge covariant) anisotropic similarity. The resulting transformation rules for the 
gauge covariant derivatives of iplx) and iplx) take the form 

Vy = RQ'^^Vip, V^V = ip VQ-^R+, (30) 

where the superscripts and ^ stand for the matrix transposition and Hermitian conjugation, 
respectively. 

2. Consider two material systems, say the dashed and the undashed ones, described by 
the thermodynamical functionals (jlj) with the densities F'{.) and F{.), respectively. We shall 
say that those systems are anisotropically similar if 

F'(T, X, ifj, tp, dtp, %p (9, A, dA) = F{T, x', tp', tp', d'tp', tp'd\ A', d'A). (31) 

Due to the gauge covariance this is equivalent to 

F'(T, X, W, V, 5) = F(T, x', ip\ Vy , 7/>V, B'). (32) 

3. To justify this conclusion let us examine the behaviour of the field equations ffT71 - [2Ul) 
under the trasformations fl25]) . Taking into account the transformation rules (1251 - [5D|) we 
obtain from eqn. (!32|) 

W~ W ^^^^ 



(9(^'V0 ^ d{ipW) 

In consequence, from the expression (120|) one obtains 



3' = Q3- (35) 
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Further, 



B' 



{det Q)-^QB 



(36) 



and 



H' 



{det Q)Q-^H. 



(37) 



The relations (133) - [37]) guarantee the covariance of the field equations (ITTl) and (ITSll with 
respect to the anisoptropic similarity transformations. 

In the layered hybrid model a straightforward calculation shows that the energy of the 
Josephson's coupling (1261) is invariant with respect to transformation 



with all the remaining quantities kept fixed. 

This transformation switches between systems with all the constants 7^ transformed 
according to the parity of q: the constants remain identical for even g's and change sign 
for odd g's; in particular = —71 and 72 = 72. At the same time the order parameter 
ipn is transformed by the sign-alternating factor; in particular the states of uniform order 
parameter are transformed into alternating ones, and vice versa. 

In consequence the energy of the Josephson's coupling is invariant with respect to trans- 
formation ( 138|) . This invariance extends to the Josephson part of the total free energy func- 
tional ( p5|) . Two hybrid systems related by the transformation ( 138|) have identical ground 
state energies, excitation energy spectra, and all the thermodynamical properties, including 
Tc- The shapes of the corresponding ground states are, however, substantially different. 

The usefulness of the above transformations stems from the fact that they allow to ob- 
tain a solution valid for one system from a solution valid for another - anisotropically similar 
- system. Analogous transformations, known as anisotropic scaling, have been defined in 
reference [5]. However, although very useful as a tool for getting approximate estimations, 
they do not imply exact invariance. In opposition to that, our anisotropic similarity trans- 
formations are exact. In consequence, transformed exact solutions for the depart system are 
always guaranteed to be also exact solutions for the target system. 
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